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It is shown that a Weakly Interacting Massive dark matter Particle (WIMP) interpretation for 
the positron excess observed in a variety of experiments, HEAT, PAMELA, and AMS-02, is highly 
constrained by the Fermi/LAT observations of dwarf galaxies. In particular, this paper has focused 
on the annihilation channels that best fit the current AMS-02 data (Boudaud et ah, 2014). The 
Fermi satellite has surveyed the 7-ray sky, and its observations of dwarf satellites are used to place 
strong bounds on the annihilation of WIMPs into a variety of channels. For the single channel case, 
we find that dark matter annihilation into {66, e + e“, p + /x“, t + t“, 4-e, or 4 -r} is ruled out as an 
explanation of the AMS positron excess (here 6 quarks are a proxy for all quarks, gauge and Higgs 
bosons). In addition, we find that the Fermi/LAT 2a upper limits, assuming the best-fit AMS-02 
branching ratios, exclude multichannel combinations into 66 and leptons. The tension between the 
results might relax if the branching ratios are allowed to deviate from their best-fit values, though 
a substantial change would be required. Of all the channels we considered, the only viable channel 
that survives the Fermi/LAT constraint and produces a good fit to the AMS-02 data is annihilation 
(via a mediator) to 4-/r, or mainly to 4-/.1 in the case of multichannel combinations. 


Weakly Interacting Massive Particles (WIMPS), such as the lightest supersymmetric particles (for reviews, see 
Refs. mm), are thought to be the best motivated dark matter (DM) candidates. The particles in consideration 
are their own antiparticles; thus, they annihilate among themselves in the early universe and naturally provide the 
correct relic density today to explain the dark matter of the universe. This same annihilation process takes place 
in the present universe wherever the DM density is sufficiently high and is the basis for DM indirect detection 
searches. Indirect detection experiments search for the annihilation products of dark matter particles, including 
electrons and/or positrons, antiprotons, photons, and neutrinos. Promising sites for the observation of dark matter 
annihilation products include the core of the Sun [3j, the Earth mm, our Galactic halo SBH31, Galactic center flOj . 
dwarf satellite galaxies mm, and from DM substructures mm- 

Over the past several years, there have been a number of experimental signals which have been interpreted as 
possible indications of DM. Confirmation that any of these observations are actually due to DM, rather than being 
a mere experimental artifact or astrophysical background, would likely require more than one experiment to provide 
complementary information. In this paper, we consider the anomalous features in the spectrum of cosmic ray positrons 
and electrons reported by AMS-02 El. PAMELA [T51 , tSj, and the Large Area Telescope of the Fermi Gamma Ray 
Space Telescope (Fermi/LAT) [20] (as well as in earlier indications from HEAT [2T1 - I23] ). The positron fraction was 
found to be a steadily increasing function of energy, above 10 GeV. This behavior is difficult to explain with standard 
astrophysical mechanisms, in which positrons are secondary particles, produced in the interactions of primary cosmic 
rays during the propagation in the interstellar medium. These observations have led to a great deal of speculation that 
DM annihilations [241143] or decays [44k47l may be responsible. However, any explanation of these positron/electron 
signals in terms of DM annihilation requires somewhat nonstandard WIMP properties. In particular, the local 
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halo density of DM within the vicinity of the Solar System is insufficient to produce these observations unless the 
annihilation cross section is considerably larger than that typically expected for a thermal relic, or the annihilation rate 
is otherwise supplemented by a large boost factor ~ 10 1 — 10 4 . Such an enhancement could arise due to astrophysics; 
for example, due to substructures in the DM distribution; yet Ref. 48] argues that the probability of such a nearby 
DM clump is < 1%. In this paper, we will instead focus on the possibility of an annihilation cross-section that is 
enhanced by the required boost factor compared to the standard thermal annihilation. Furthermore, only a handful 
of possibilities can explain the spectral shape reported by AMS-02 and PAMELA, as well as avoid overproducing 
antiprotons (in excess of what is observed) l4Tjl - f4?l . The most well-studied approach to satisfying all these constraints 
has been leptophilic DM, i.e. the DM annihilations proceed largely to leptons or t~t~ ), which do not produce 

any antiprotons. Alternatively, other annihilations channels—quarks, vector and Higgs bosons—are allowed if the 
DM particle is heavier than ~10 TeV, producing antiprotons at higher energies than those probed by AMS-02 data 

US- 

Other explanations of the Cosmic Ray Positron Excess (henceforth CRPE) have also been proposed. The most 
plausible is that it is due to pulsars [53ti56j . 

Positrons can ultimately be produced by DM annihilating into a variety of channels: single channels (leptons, 
quarks, gauge bosons, Higgs, or four leptons) or a mix of these channels. Prior to 2014, when the most recent AMS-02 
results became available, all annihilation channels could provide explanations for the positron excess (see Table 2 
in Ref. [57]). The most popularly studied cases were WIMP masses ~200 GeV, which matched the data only for 
leptophilic channels (DM annihilates only into leptons or t + t“), though quark and gauge boson channels well 

fit the data at higher WIMP masses, with best-fit masses as high as ~50-200 TeV (depending on the channel). 

The recent AMS-02 data release m has greatly improved our understanding of the positron excess. As stressed 
by Ref. [57;, two major improvements have emerged. First the new data are far more accurate and extend out to 
500 GeV, much higher energies than previously explored. Improved accuracy in the positron spectrum leads to stronger 
constraints on any model for the origin of the positrons. Second, AMS-02 has measured directly the total electron and 
positron flux, the denominator in the positron fraction, reducing systematic errors and leading to differences in best-fit 
regions for the DM mass and cross section. Ref. m further stressed the importance of addressing uncertainties in 
the cosmic ray (CR) propagation model, as these models are of critical importance in assessing the true nature of the 
positron excess. 

Boudaud et al. m have performed an analysis with the new AMS-02 data, exploring the DM models that could 
explain the CRPE. Compared to earlier data and analyses, the recent AMS-02 data do indeed lead to much stronger 
constraints on suitable DM candidates that can explain the e + excess and the best-fit masses have shifted by a fair 
amount. Primarily these authors used the benchmark set of CR propagation model parameters known as MED in 
obtaining their results; this is the model that best fits the B/C ratio in the cosmic rays. With the recent, more 
accurate data and using MED, these authors found that many of the earlier allowed models have been ruled out (see 
Table 1 in Ref. [57]). In terms of single channel models, leptophillic DM is no longer viable, whereas single channel 
annihilation into quarks and gauge bosons provide an excellent fit to the CRPE. They also considered the case of 
DM annihilation into a mediator particle </>, where \X with </> —> ££, thus yielding four leptons per annihilation 

(“4-lepton” channel). Again, single channel annihilation into four leptons does not fit the AMS-02 data (assuming 
MED propagation). Alternatively, annihilation into a combination of channels can also provide an excellent fit, in 
particular annihilation into an admixture of leptons and bb pairs. Finally, a combination of four lepton channels 
(arising from a mediator field </>), specifically the four-tau (75%) and four-electron (25%) channel, turns out to provide 
a good fit to the AMS-02 data for a DM mass between 0.5 and 1 TeV. The results reprised in this paragraph all refer 
to the MED CR propagation model. 

Far larger than the statistical errors in the AMS-02 data are the systematic errors associated with the CR propaga¬ 
tion. Thus Boudaud et al. m looked at a set of 1623 different combinations of the CR transportation parameters—all 
consistent with observed boron-to-carbon ratios—which bracket the systematic uncertainty in the propagation model 
(see also Ref. [SB] for an analysis of the new AMS-02 data, with a different treatment of the propagation). By including 
these propagation model parameters into the fits to the AMS-02 data, specifically by finding the DM parameters that 
fit the data for any of the 1623 propagation parameter combinations, the allowed DM parameter space increases. In 
terms of the single channel models, the 2-r, 4-/i, and 4-r cases can now provide excellent fits to the positron excess 
in addition to the cases mentioned above. 

DM annihilations that produce e + , either directly or through decays and showering of the primary annihilation 
products, will invariably also produce y-rays. Thus we can constrain the above DM models by comparing with y-ray 
observations. The best places to look are regions with a large abundance of DM: the Galactic Center, clusters, and 
dwarf galaxies. We will focus on dwarf galaxies for the remainder of the paper. The dwarf spheroidal galaxies inside 
the Milky Way are some of the most dark matter dominated objects known, with mass-to-light ratios as high as ~ 
1000. Because they are so dark matter rich and nearby, they are exceptionally good places to indirectly detect DM via 
y-rays produced in its annihilation. For previous y-ray data that could be used to constrain a DM interpretation of a 
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positron excess, see the previous Fermi/LAT combined analysis of dwarf galaxies [55], VERITAS [60] . and MAGIC EP 
observations of Segue 1, as well as results of the H.E.S.S. collaboration [62] on Sagittarius and other dwarf galaxies. 

Fermi/LAT surveys the y-ray sky and, specifically, has looked for y-rays from 25 dwarf galaxies, detecting no 
significant excess [55] . This lack of signal is used to place y-ray flux upper limits for energies between 500 MeV and 
500 GeV. These bounds are then used to constrain the dark matter annihilation for a broad range of particle masses 
and annihilation channels. 

The differential y-ray flux from DM annihilation in a dwarf galaxy can be written as the product of two 
components, a factor ]^ pp that encodes all the particle physics and the so-called J-factor that contains the astrophysics, 
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Here, ( av) is the DM annihilation cross section, m x is the DM mass, and jA is the photon spectrum from the DM 
annihilation, which depends on the DM mass m x and the annihilation channel. The J-factor (the term in the second 
set of parentheses) integrates the square of the DM density p x along the line of sight and over a solid angle Af2. The 
J-factor can be estimated from stellar kinematics as stars act as tracers of the gravitational potential, allowing the 
DM distribution to be inferred. 

We use Fermi/LAT dwarf y-ray results to constrain DM models using the data and likelihood technique described 
by the Fermi/LAT collaboration in Ref. [55], to be briefly reviewed here. The primary quantity used in the likelihood 
analysis is the energy flux 
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for each dwarf (indexed by k) and energy bin (indexed by j). Here, ^ = g§ pp x Jk is the differential flux for a dwarf 
with J-factor Jk- For each dwarf and energy bin, Fermi/LAT provides a likelihood Ckj in s k,j- The likelihood function 
accounts for instrument performance, the observed counts, exposure, and background fluxes. For a given annihilation 
channel, the energy flux is dependent only on the theoretical parameters m x , {av ), and Jk] he. Sk,j = Skj{m x , {av), Jk)- 
Accounting for observational constraints on the J-factor, the likelihood for a given dwarf Cu is 


Ck{m x ,{av),J k ) = CJ\[(Jk\Jk,ak)~WC-k,j(s k ,j{m x ,{av), Jk)) , 
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(3) 


where CAf represents a log-normal distribution and Jk & ak are the parameters describing that distribution, derived 
from the stellar kinematics in the dwarf. The combined likelihood for multiple dwarfs is 


£(m x ,(av),J) =Y[Ck(m x ,(crv),Jk), (4) 

k 

where J represents the set of J-factors { Jk}, one for each dwarf. We use in our analysis the same 15 non-overlapping 
dwarf galaxies with J-factor estimates that are used by Fermi/LAT in their analysis. We use the J-factor estimates 
based upon a Navarro-Frank-White (NFW) density profile [551, though a Burkert profile [53] would not significantly 
affect our results. 1 

Fermi/LAT constraints in {av) at a given m x are determined using a delta-log-likelihood approach treating the 
J-factors as nuisance parameters. The delta-log-likelihood AC is given by 

Aln£(m x , {av)) = In C(m x , {av),J) — ln£(m x , {av), J) , (5) 

where {av) & J are the values of {av) & J that jointly maximize the likelihood at the given m x and J = J(m x , {av)) 
are the J-factors that maximize the likelihood for the given m x and {av). The ID confidence intervals in {av) at the 
na confidence level are determined by identifying the range of {av) such that 

A In C{m x , {av)) < n 2 /2 . (6) 


1 References [591 1651 have shown that the the integrated J-factor within 0.5 degrees is fairly insensitive to the choice of dark matter 
density profile so long as the central value of the slope is less than 1.2. 
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FIG. 1: Constraints on the dark matter annihilation cross-section ( av) and mass m x for annihilations into bb (left) and 
four r’s via a mediator <j> (right). The best-fit AMS-02 parameters, as derived by Ref. [57] for the MED propagation model 
parameters, are shown in red, while the Fermi/LAT upper bounds are shown by the blue curves. Fermi/LAT dwarf constraints 
are generated using the procedure described in Ref. [59] • Constraints are shown at the 2a confidence level. For the bb case, the 
AMS-02 best-fit points for a selection of other cosmic ray (CR) propagation model parameters are shown as gray dots (also 
taken from Ref. 1571 1. 


We will generally show the upper limit of the 2 er confidence intervals (95.4% confidence level). 

We consider the case of DM annihilating directly to bb and leptons ££, as well as the “4-lepton” channels (via a 
mediator </>). We use the bb case as a proxy for all other quarks, gauge bosons, and the Higgs boson. The spectra for 
the u, d , c, s, and t quarks; the W, Z , and g gauge bosons; and the H boson are all similar in shape to the b spectrum 
and within ^50% of the amplitude. Such differences of at most a factor of two in amplitude will not prove to be 
significant; thus the b spectrum should be reasonably representative of these other cases . 2 The spectra from leptons, 
on the other hand, depend on the flavor, so we consider e, p,, and r separately. The annihilation spectra ^ for these 
different channels are derived using PYTHIA8 El EH- For the first set of cases, annihilation directly to quarks or 
leptons, we include final state radiation (FSR) in the PYTHIA simulations as these are the primary source of photons 
for the e + e _ and p + p~ channels; the FSR provides only a minor contribution to the bb and t + t~ channels. We do 
not include FSR for the mediator case due to the complexity, though that is again the primary source of photons for 
the 4-e and 4-/x channels. 

The AMS-02 2er confidence regions in (av) vs. m x for the bb and 4-r channels are shown in red in the left and 
right panels, respectively, of Figure |TJ taken from Ref. m- The Fermi/LAT 2er upper limits in (av) are shown by 
the blue curves in the figure. For both channels, the AMS-02 regions are strongly excluded by the Fermi/LAT dwarfs 
data. In Figure [l] the AMS-02 2cr confidence regions were determined assuming the MED propagation parameters, 
the set of five astrophysical parameters fixed to best fit the measured B/C ratio EH] Hi • As the choice of propagation 
parameters has an impact on the fit to the AMS-02 data, we show also the best-fit points found by Ref. [57] for a 
selection of other CR propagation model parameters as discussed above. These best-fit points, shown as gray dots for 
the bb case, are taken from Figure 14 of Ref. [57] (parameters with p > 0.0455 only). Though varying the propagation 
parameters broadens the region of parameter space consistent with the AMS-02 data, the DM interpretation of the 
positron spectrum remains strongly in conflict with the Fermi/LAT dwarf 7 -ray observations. 

We wish to comment on the differences between regions compatible with AMS-02 assuming the MED propagation 


2 Each of the AMS-02 and Fermi constraints vary by less than a factor of two from the bb case when looking at other quarks, but bb is 
ruled out by an order of magnitude, so a factor of two is not enough to evade the Fermi/LAT constraints. 





5 



m x [TeV] 


FIG. 2: AMS-02 best-fit parameters and Fermi/LAT constraints on (av) and m x for the bb and li channels. Best-fit AMS- 
02 values for single annihilation channels are shown as red circles (taken from Ref. |57| 1. Empty circles are for the MED 
propagation parameters and filled circles are for the best-fit propagation parameters. Fermi/LAT constraints for each of these 
single channels are shown as solid blue lines. The red and blue dashed curves represent respectively the AMS-02 best-fit and 
Fermi/LAT constraint on (av) at each mass for the mixed bb+£l channels, assuming a common (uniform) branching ratio (BR) 
into each of the three leptons; Fermi/LAT constraints are generated assuming the AMS-02 best-fit BRs at each mass. The 
dotted curves are the same, but relaxing the assumption of uniform lepton BRs. We note that the MED propagation parameters 
are used for the mixed channels. While we do not believe that changing the CR propagation parameters will significantly affect 
our results, proving this assertion is beyond the scope of this paper. 


parameters and the regions compatible when the propagation parameters are allowed to vary over reasonable values, 
as represented by the 1623 sets of CR propagation parameters. In Figure [T] one can see that all the points, for both 
MED and other CR propagation model parameters, are in an elongated region roughly parallel to the Fermi/LAT 
bounds. Thus changing the CR propagation parameters does not alleviate the tension between Fermi/LAT and a 
DM interpretation of AMS-02. We believe that this statement will be generally true, not just for the bb channel but 
for all channels (single and mixed). 3 Thus we believe the following: if a DM annihilation channel that fits AMS-02 
with MED propagation is ruled out by Fermi/LAT, then we believe that the same channel is also likely to be ruled 
out if other propagation parameters are used. In other words, we do not believe that using different reasonable CR 
propagation parameters will change our results. However, we cannot prove this assertion without a detailed reanalysis 
of the AMS-02 data, beyond the scope of the current paper. 

Figure [2] shows the AMS-02 best fits and Fermi/LAT dwarf constraints on (av) and m x for annihilations into single 
channel be or bb, as well as mixed cases £1 + bb. Note that Figure [2] is over a broader mass range than Figure [l] 
First let us discuss the single channel cases. The Fermi/LAT upper bounds on the single channel cases are shown 
as solid blue lines, from bottom to top: bb, t + t~ , y + yT, e + e~. The e and y, channels have weaker limits as photon 
production is suppressed in these cases, here coming only from the FSR. 4 The b and r channels produce photons 
through unsuppressed shower/decay processes; the presence of FSR has little impact on the constraints for these two 
channels. The best-fit AMS-02 points for these annihilation channels are shown by the red circles. Empty circles 
are for the MED propagation parameters and filled circles are for the best-fit propagation parameters [57] : the DM 


3 Since the annihilation rate scales as (<jv)/m x , moving along the elongation line corresponds to roughly a fixed number of e+, as is 
required to match the AMS-02 data. There is a slight effect on the fit due to the m x dependence of the spectrum, but the general trend 
still holds. 

4 The muon produces a photon in 0(1%) of its decays, though this process is not accounted for in PYTHIA. 
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FIG. 3: Same as Figure [5] but for the mediated 4-lepton channels. The solid blue line is the Fermi/LAT bound for the 4-r 
case, which excludes the best-fit AMS-02 values. The 4-e and 4 -fi Fermi/LAT bounds have not been calculated here, but are 
likely to be weaker than their corresponding best-fit AMS-02 parameters (see the text). Whereas the 4-e case does not provide 
a good fit to the AMS-02 data, the 4-/r case does and thus survives as a viable alternative. 


parameters do not differ much between the two cases for any of the channels. The AMS-02 b and r best-fit points 
are excluded by the Fermi/LAT constraints, while the e and /i points are not. The two lighter lepton channels fail 
to be excluded by Fermi/LAT because they provide positrons in abundance for the AMS-02 signal, while photon 
production is suppressed, leading to little expected sensitivity to these channels via y-ray searches. However, both 
the e and /r annihilation channels, while capable of producing substantial numbers of positrons, are simply a poor 
fit to the AMS-02 spectrum. Thus, none of these four channels can provide a reasonable fit to both the AMS-02 
and Fermi/LAT results and no single annihilation channel into quarks, leptons, or gauge bosons can simultaneously 
explain the AMS-02 data while remaining in agreement with Fermi-LAT bounds. We note that both the MED and 
the best-fit CR propagation parameters have been used in studying these single channel cases. 

Reference m also considered the mixed-channel case, where the DM annihilates into some combination of bb 
and the three leptons. Here the MED propagation model is used for the mixed channels. Two possibilities were 
considered: one where all three leptons were assumed to have a common (uniform) branching ratio (BR) and one 
where this assumption is relaxed. The AMS-02 best-fit {av) as a function of mass are shown as red dashed and dotted 
curves for the uniform and non-uniform lepton cases, respectively (taken from Figures 7 & 5 in their paper). Their 
results were presented only down to DM masses of 0.6 TeV, hence the termination of the curves at that mass. For the 
uniform case, the leading annihilation channel is always to bb, while for the non-uniform case, t + t~ dominates below 
20 GeV and bb dominates above. The corresponding Fermi/LAT 2 a upper limits, assuming the best-fit BRs, are 
shown in blue dashed and dotted curves. The AMS-02 best-fit {av) and BRs are strongly excluded by the Fermi/LAT 
results. This does not rigorously imply the AMS-02 and Fermi/LAT results are in strong conflict for all multi-channel 
cases, as the Fermi/LAT constraints will vary if the BRs are allowed to deviate from their best-fit values. However, 
bringing the two results into compatibility will require the bb channel to be heavily suppressed, as well as the r + r _ 
channel at the lighter end of the mass range, which is quite different from the best-fit case, where bb dominates at 
higher masses and bb + t + t~ accounts for > 50% of the annihilations at lower masses. 

Figure [3] examines the mediated 4-lepton case with \X -► </ >( / ) and <t> The 4-r Fermi/LAT upper limit is 

shown by the solid blue curve. Due to the complexity of implementing their non-leading-order photon production 
mechanisms, we have foregone the calculation of limits for the 4-e and 4 -/i cases. However, as with the 2-lepton cases 
shown in Figure [2} these e and /r limits can be expected to be 1-2 orders of magnitude weaker than the r limit. The 
AMS-02 best-fit points are again shown by red circles, determined using both MED propagation (empty circles) and 
the best-fit propagation parameters (solid circles) [57]. The r point is excluded by the Fermi/LAT results. Without 
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a calculation of the e and fi Fermi/LAT constraints, compatibility of the two experimental results cannot be checked, 
though the Fermi/LAT constraints are almost certainly too weak to exclude these two points. The e channel is a poor 
fit to the AMS-02 data and thus of little interest. However, the channel both provides a good fit to the AMS-02 
data and escapes Fermi/LAT constraints, thus remaining a viable model. 

Reference [571 considered the mixed 4-lepton case, with the best-fit (av) shown by the dotted red curve in Figure [3j 
using the MED propagation model. The best-fit BRs are annihilation exclusively to taus for DM masses above 2 TeV, 
and a mix of electrons and taus below that, though still dominated by the tau channel for the masses shown. The 
Fermi/LAT constraint for these best-fit BRs is shown by the dotted blue curve, which becomes identical to the single¬ 
channel 4-r constraint above 2 TeV. The AMS-02 best-fit (av) and BRs are incompatible with the Fermi/LAT data 
except at very high masses ( m x > 40 TeV). The analysis shown in Figure [3] includes the case featured by Ref. El] as 
their favored scenario: a combination between the four-tau (75%) and four-electron (25%) channels for a DM mass 
between 0.5 and 1 TeV. This case is ruled out by Fermi/LAT as an explanation of AMS-02 data (for m x < 40 TeV). 
The caveats discussed for the bb+££ case previously apply here: allowing the BRs to vary from their best-fit values 
will change the Fermi/LAT constraints and could potentially bring these experimental results into line for this model. 
The t channel would need to be strongly suppressed for that to occur, which implies BRs far different than the best-fit 
case. An example where this happens is the predominantly 4-/r case which remains a viable model. 

In addition to bounds from 7 -rays and antiprotons already discussed, there are other observational constraints on 
a DM interpretation of the e + excess seen in AMS-02. DM annihilation products also include neutrinos. References 
[70HZ2] pointed out that a leptophilic explanation of a positron excess should also produce large numbers of neutrinos 
detectable in the IceCube neutrino observatory. Currently IceCube upper limits are in tension with the DM bestfit to 
the positron anomaly, assuming MED propagation parameters, when DM annihilates into W + W~ 73). In addition, 
measurements of the CMB temperature and polarization provide constraints on the annihilation cross section of DM 
EMS]. In the near future, papers from the Planck Satellite with extremely strong bounds are expected and could 
rule out the DM scenario completely. 

In summary, we have used the Fermi/LAT dwarf galaxy data in order to constrain dark matter as an explanation of 
the positron excess seen in HEAT, PAMELA, and AMS-02. In particular, this paper has focused on the annihilation 
channels that best fit the current AMS-02 data E3- We first considered the single channel case of DM annihilating 
directly to bb and leptons ££, as well as to 4-leptons (via a mediator </>); we then considered the multi-channel case 
where annihilation proceeds through a combination of channels. We used bb as our proxy for all other quarks, gauge 
bosons, and Higgs bosons, since the spectra and amplitudes are similar, but considered each lepton flavor separately. 
For the single channel case, we found that dark matter annihilation into {bb, e + e _ , t + t~, 4-e, or 4-r} cannot 

both provide a good fit to AMS-02 and avoid the 2er upper limit from Fermi/LAT. The AMS-02 best-fit BRs and (av) 
used for the analysis of this paper assume either the MED propagation model defined in Refs. [ 68 ] [69] or the best-fit 
propagation parameters E3- Multi-channel annihilations are also highly constrained by Fermi/LAT’s measurement 
of 7 -ray flux. We use the MED propagation model parameters for the mixed channel cases, but doubt that other 
reasonable choices of CR propagation parameters would change our results. Specifically, we find that the Fermi/LAT 
2a upper limits, assuming the best-fit AMS-02 BRs, exclude the annihilation into a combination of bb and the three 
leptons for DM masses 600 GeV < m x < 100 TeV. In addition, the Fermi/LAT upper limit is incompatible with the 
AMS-02 best fit (av) and BRs for annihilation into a mix of the mediator driven 4-lepton channels except for very 
high masses (m x > 40 TeV). However, this does not rigorously imply that the results from AMS-02 and Fermi/LAT 
are in strong conflict for all multi-channel cases, as the Fermi/LAT constraints will vary if the BRs are allowed to 
deviate from their best fit AMS-02 values. In order to reconcile both experiments, the BRs considered would have to 
deviate significantly from their best-fit values. We find that the dark matter annhiliation into 4 -[j, provides a good fit 
to the AMS-02 data and escapes Fermi/LAT constraints. The reason for the 4-/r channel escaping Fermi/LAT upper 
limits is that it provides positrons in abundance for the AMS-02 signal, while photon production is suppressed. Hence, 
for the best-fit BR to the AMS-02 data, we find only one viable DM annihilation channel that survives Fermi/LAT 
constraints and provides a good fit to the current AMS-02 data: the 4-/x channel. 

We briefly mention ways around these conclusions. First, for the mixed channel cases we have only studied the best 
fits to AMS-02 data (found in Ref. [57]). It is possible that there are other BRs that still provide reasonably good fits 
(though not the best fits) to the AMS-02 data that are not ruled out by Fermi/LAT bounds from dwarf spheroidals. 
A joint statistical analysis of the AMS-02 and Fermi/LAT would be required to check for other alternatives. Second, 
although it is extremely unlikely that the boost factor required by the AMS-02 data is due to a nearby clump of DM 
that is 100-1000 times as dense as its surroundings, perhaps part of the boost factor is due to a clump of, say, a 
factor of 10 (again, unlikely). In that case the required annihilation cross-section to explain the AMS-02 data could be 
lower, and more channels would remain viable. Third, it is possible that part of the AMS-02 signal is due to pulsars, 
and part due to DM annihilation. Again, more DM annihilation channels might then still remain compatible with 
Fermi/LAT bounds from dwarf spheroidals. These latter two caveats would also apply to bounds on these scenarios 
from the CMB including those expected from upcoming Planck data. 
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